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We study the behavior of the pion dispersion relation in a pion medium at finite density and 
temperature. We introduce a pion chemical potential to describe the finite pion number density and 
argue that such description is valid during the hadronic phase of a relativistic heavy-ion collision 
between chemical and thermal freeze-out. We make use of an effective Lagrangian that explicitly 
respects chiral symmetry through the enforcement of the chiral Ward identities. The pion disper- 
sion relation is computed through the computation of the pion self-energy in a non-perturbative 
fashion by giving an approximate solution to the Schwinger-Dyson equation for this self-energy. 
The dispersion relation is described in terms of a density and temperature dependent mass and an 
index of refraction which is also temperature, density as well as momentum dependent. The index 
of refraction is larger than unity for all values of the momentum for finite p and T. We conclude 
by exploring some of the possible consequences for the propagation of pions through the boundary 
between the medium and vacuum. 
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I. INTRODUCTION 

The phase structure of hadronic matter at high den- 
sity and temperature has become a subject of increasing 
attention in recent years, both from the theoretical and 
experimental points of view. This attention is driven 
mainly by the possibility to produce a locally thermal- 
ized, deconfined state, whose degrees of freedom arc the 
quarks and gluons of quantum chromodynamics (QCD), 
the so called quark-gluon plasma (QGP), in high-energy 
heavy- ion collisions Q. 

If the QGP is produced in these kind of reactions, the 
prevailing view portraits the evolution of such a system 
traversing a series of stages, the last of which consists of 
a large amount of hadrons strongly interacting in a finite 
volume, until a final freeze-out. 

For not too high temperatures, hadronic matter con- 
sists mainly of pions, therefore, the study of the prop- 
agation properties of pions within the above described 
environment represents an important ingredient for the 
understanding of the properties of the hadronic system at 
and just before freeze-out. For instance, it has been spec- 
ulated that the pion interactions within a dense hadronic 
environment can give rise to interesting collective surface 
phenomena ||. Moreover, it is also well known that the 
pion group velocity is an important piece of information 
necessary to properly account for the dilepton spectrum 
coming from the hadronic phase |3|. 

The hadronic degrees of freedom are customarily ac- 
counted for by means of effective chiral theories that in- 
corporate the Goldstone boson nature of pions. One of 
such theories is chiral perturbation theory (ChPT) which 



has been employed to show the well known result that 
at leading perturbative order and at low momentum, 
the modification of the pion dispersion curve in a pion 
medium at finite temperature is just a constant, temper- 
ature dependent, increase of the pion mass ChPT 
has also been used in a two-loop computation of the pion 
self-energy [|| and decay constant ||. A striking result 
obtained from such computations is that at second order, 
the shift in the temperature dependence of the pion mass 
is opposite in sign and about three times larger in mag- 
nitude than the first order shift, already at temperatures 
close to 150 MeV. This result might signal either the 
breakdown of the perturbative expansion at these tem- 
peratures or the need to perform such calculations using 
schemes other than perturbation theory. 

However, a missing ingredient in the calculations of the 
pion dispersion curve is the treatment of the medium's 
finite density. The conceptual difficulty is related to the 
fact that, though it is possible to assume that the system 
is in (at least local) thermal equilibrium, strictly speak- 
ing the only conserved charge that can be associated to 
the pion system is the electric charge and thus, for an 
electrically neutral pion system the corresponding chem- 
ical potential vanishes. The behavior of the pion mass 
in the presence of an isospin chemical potential has been 
recently studied in Ref. j7j. But in order to describe a 
situation in which the number of pions in thermal equi- 
librium is finite, we need to consider a chemical potential 
associated with the pion number, instead of its charge. 

Recall that the pion number is not a conserved quan- 
tity due to either strong, weak or electromagnetic pro- 
cesses. Nevertheless, the characteristic time for electro- 
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magnetic and weak pion number-changing processes, is 
very large compared to the lifetime of the system created 
in relativistic heavy-ion collisions and therefore, these 
processes are of no relevance for the propagation proper- 
ties of pions within the hadronic phase of the collision. 
As for the case of strong processes, it is by now accepted 
that they drive pion number toward chemical freeze-out 
at a temperature considerably higher than the thermal 
freeze-out temperature and therefore, that from chemi- 
cal to thermal freeze-out, the pion system evolves with 
the pion abundance held fixed [|[ ||. Under these cir- 
cumstances, it is possible to ascribe to the pion density a 
chemical potential and consider the pion number as con- 
served [[l0[ In this context, the role of a finite pion 
chemical potential into a hadronic equation of state has 
been recently investigated in Ref. Jl2|. 

Furthermore, another important ingredient in the anal- 
ysis is the well know fact that in finite temperature field 
theories with either massless degrees of freedom or that 
exhibit spontaneous symmetry breaking Jl3| , the pertur- 
bative expansion breaks down and thus the necessity to 
implement resummation techniques. 

In this paper we explore the effects introduced by a 
finite pion density on the dispersion curve of pions at fi- 
nite temperature. Starting from the linear sigma model, 
we use an effective Lagrangian jl4| obtained by integrat- 
ing out the heavy sigma modes and compute, in a non 
perturbative fashion, the pion self-energy. We find that 
the dispersion curve is modified with respect to the vac- 
uum in a way described by the introduction of an index 
of refraction larger than one, in addition to the thermal 
and density increase of the pion mass. We discuss possi- 
ble implications of such behavior of the pion dispersion 
curve. 

The work is organized as follows: In section |n[ we 
recall how for temperatures and momenta smaller than 
the mass of the sigma, it is possible to integrate out the 
heavy sigma modes to construct an effective Lagrangian 
that reproduces the two- loop pion self-energy at finite 



temperature. In section [II, we use this effective La- 



grangian to compute non-perturbatively the pion self- 
energy and from it, the pion dispersion curve at finite 
density and temperature. Finally, in section IV, we sum- 



marize and discuss our results, emphasizing some of the 
possible consequences of the behavior of the dispersion 
curve thus found for the propagation pions that approach 
the boundary between the hadronic medium and vacuum 
in a relativistic heavy-ion collision. 



II. EFFECTIVE LAGRANGIAN 

The Lagrangian for the linear sigma model, including 
only meson degrees of freedom and with an explicit chiral 
symmetry breaking term, can be written as [p_5| 
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FIG. 1: Diagramatic representation of the self-consistent pion 
self-energy. The double lines represent the full pion propaga- 
tor. 
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where tt and a are the pion and sigma fields, respectively, 
and the coupling A 2 is given by 
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From the above Lagrangian one obtains the Green's 
functions and Feynman rules to be used in perturbative 
calculations, in the usual manner. For instance, the bare 
pion and sigma propagators A w (P) , A CT (Q) and the bare 
one-sigma two-pion and four-pion vertices r^ 2 , are 
given by (hereafter, capital Roman letters are used to 
denote four momenta) 
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These Green's functions are sufficient to obtain the mod- 
ification to the pion propagator, both at zero and finite 
temperature, at any given perturbative order. 

When interested in a given approximation scheme to 
build such Green's functions, it is possible to exploit 
the relations that chiral symmetry imposes among them. 
These relations, better known as chiral Ward identities 
(ChWI), are a direct consequence of the fact that the di- 
vergence of the axial current may be used as an interpo- 
lating field for the pion M|. Thus, one can construct the 
modification to one of the above Green's functions at a 
given perturbative order and for a given approximation. 
In order to make sure that the approximation respects 
chiral symmetry, one needs to check that the modifica- 
tion of other Green's functions respect the corresponding 
ChWI. For example, two of the ChWI satisfied -order by 
order in perturbation theory- by the functions A 7r (P), 
A a (Q), T% and T& kl are 

/ 7r r^( ; o,Pi,p 2 ,p 3 ) = r«(p i; p 2 ,p 3 )^ 

+ T l ^(P 2] P 3 ,P 1 )8 ik 

+ T{ k (P 3 ;P 1 ,P 2 )S il 

/,r«(Q;0,P) = [A- 1 (g)-A- 1 (P)]^(4) 

where momentum conservation at the vertices is implied, 
that is Pi + P 2 + P 3 = and Q + P = 0. The notation 
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for the functional dependence of the vertices in Eqs. (Q) 
is such that the variables before and after the semicolon 
refer to the four-momenta of the sigma and pion fields, 
respectivel y fl5|l . 

In Refs. |14| it has been shown that in the kinematical 
regime where the pion momentum, the pion mass and 
the temperature are small compared to the sigma mass, 
the effective one-loop sigma propagator and one-sigma 
two-pion and four-pion vertices are given by 
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x [1 - 3X 2 l t (P 1 +P 2 )] S lj 5 kl 
+ [l-3A 2 X*(P 1 +P 3 )] S ik S jl 
+ [l-3A 2 J*(Pi + P 4 )] 5 il S jk }, 

(7) 

where in the imaginary-time formalism of thermal field 
theory (TFT), the function X* is obtained as the time- 
ordered analytical continuation to real energies of the 
function X defined by 



K 2 + ml {K - Q) 2 + ml 
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Here Q = (w,q), K = (w n ,k), Q 2 



2 + k 2 with lj — 2mnT and i 



2nirT (to, n integers) 
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being discrete boson frequencies, T is the temperature 
and q = q|, k = |k|. 

It is easy to check that Eqs. (||), (|^) and ([?]) satisfy the 
Ward identities in Eq. (|J), this ensures that the approx- 
imation scheme adopted respects chiral symmetry. 

By using the above effective vertices and propagator, 
it is possible to construct the two-loop modification to 
the pion self-energy in the same kinematical regime with 
the result @ 

n 2 (P) = (2|) T E/ ( 2^)3 K 2 + m 2 

x {s - {^j [92*(0) + 6J f (P + *)]}• (9) 

Equation (g) reproduces the leading order result obtained 
from ChPT Q. Furthermore, we observe that Eq. (||) 
can be formally obtained by means of the effective La- 
grangian 
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where a = 6(m 2 /2f 2 ) and the factor 6 comes from con- 
sidering the interaction of like-isospin pions in the vertex 

\2Jtv J 

Eqs. (g) and © mean that in the kinematical regime 
where the sigma mass is large compared to the pion mass, 
the momentum and the temperature, the linear sigma 
model Lagrangian reduces to a 4 Lagrangian for effec- 
tive like-isospin pions with an effective coupling given 
by a. By restricting ourselves to the above kinematical 
regime, we will proceed on working with the Lagrangian 
given by Eq. ( |l0| ) . 



III. NON-PERTURBATIVE PION 
SELF-ENERGY 

It is well know that quantum field theories at finite 
temperature present certain subtleties such as the break- 
down of the perturbative expansion |]l6| . This break- 
down becomes manifest in two important cases: the ap- 
pearance of infrared divergences in theories with massless 
degrees of freedom, and the compensation of powers of 
the coupling constant with powers of T for large temper- 
atures. In both situations, the resummation of certain 
classes of diagrams represents an important improvement 
for the study of the physical properties of such theories. 
Even for cases where neither there is a massless degree 
of freedom, nor the temperature is extremely large, it is 
important to consider a resummation scheme, particu- 
larly for the case of theories with spontaneous symmetry 
breaking near critical behavior. 

In order to consider the above mentioned general situa- 
tion and with the purpose of studying the pion dispersion 
curve at finite density and temperature, let us formally 
consider the Schwinger-Dyson equation for the pion self- 
energy depicted in Fig. 1, whose analytic expression is 
given by 
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where for internal lines, we make the substitution iu> — > 
iuj + fi. Notice that since the interaction Lagrangian con- 
tains only a quartic term, there is no need to dress the 
vertices in the above equation. 

Equation (TL2) represents an integral equation for the 
function n(Pj7 which, needless to say, is very difficult 
to be solved exactly. In order to find an approximate 
solution let us write 



n(p) ~n + n(p) 



(13) 
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and consider II(P) <C H). As we will see, such as- 
sumption is justified given that in our approximation, 
II ~ 0(a 2 ). Keeping only the lowest order contribution 
in n, Eq. ( |l2] ) becomes 



n(P) = n + n(P) 



d 3 k 
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which in turn, serves as the definition of the functions IIo 
and n(P), given by 
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Equation (|l5|) represents a self consistent equation for the 
(momentum independent) constant IIo. This is the well 
known resummation for the superdaisy diagrams which 
constitute the dominant contribution in the large- N ex- 
pansion |ll| of the Lagrangian in Eq. (|To|) . On the other 
hand, Eq. (|l^) represents a first approximation to the 
momentum dependent piece of II(P). 

The solution to Eq. ( |l5| ) is given by the transcendental 
equation 
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Figure 2 shows the behavior of the quantity \J ml + IIo 
as a function of (a) the temperature T for different values 
of the chemical potential p and (b) as a function of \i for 
different values of T. In both cases, -J m\ + Hq grows 
monotonically with both T and fi. 

The calculation of the function IT(P) involves that of 
the function 



<s(p)^t£ 



d 3 k X(K;^/ml + U ) 
(2tt) 3 (K - P) 2 + ml + lie 



(18) 



where T(K; y/ml + n ) is the function defined in Eq. (g), 
with the replacement m^ — > \J ml + Hq . The calculation 
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FIG. 2: VmJ + n as a function of (a) the temperature T for 
different values of the chemical potential ranging from n — 
to n = 130 MeV and as a function of (b) the chemical potential 
for different values of the temperature ranging from T = 50 
to T = 150 MeV. 



of the dispersion relation requires knowledge of the real 
part of the retarded version of the function 5(P), which, 
after analytical continuation to real frequencies, is given 

by @ 
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where V stands for the principal part of the integral, e is 
the sign function and 
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is the Bose-Einstein distribution. 

From Eq. ([l9]) we see that we require only knowledge 
of the imaginary part of T l (ko, k). It is a lengthy, though 
straightforward exercise to show that 



ImJ t (fc ,A;): 



1 

llfrr 



e[^ 2 _4(m2 +n )] 



5 




50 75 100 

p (MeV) 

FIG. 3: p as a function of p for /j, = 0, 120, 130 MeV and 
T = 120 MeV. 
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where i-C 2 = fc 2 — fc 2 , is the step function and the 
functions a and u>± are 



a(^ 2 ) = \ I 
a; ± (fe , fc) 
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It can also be checked that for fj, — 0, Eqs. ( |21 
and ( p2[ ) reduce to the corresponding expressions found 
in Refs. 0. 

Notice that Eq. (|l^) contains temperature-dependent 
infinities coming from the terms involving the function 
Z*, as well as vacuum infinities. This is an usual feature 
of multi-loop calculations at finite temperature where one 
always encounters temperature-dependent infinities in in- 
tegrals involving only the bare terms of the original La- 
grangian. However, as it turns out, these infinities are 
exactly canceled by the contribution from the integrals 
computed by using the counterterms that are necessary 
to introduce at the previous order of the loop expan- 
sion to carry the (vacuum) renormalization. This renor- 
malization method has been developed in Refs. for 
the resummation method employed here and has been 
given the name Modified Self- Consistent Resummation. 
For the purposes of this work, we concentrate on the fi- 
nite, temperature-dependent terms and refer the reader 
to the above cited works for a detailed treatment of the 
renormalization procedure. 

According to Eqs. and 
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FIG. 4: Fit of the dispersion curve to the functional form po = 
y/n- 1 ^,^ 2 + M 2 (T, fi) for (j, = 0, 120, 130 MeV and T = 
120 MeV. The dots are the computed values of the dispersion 
relation. 
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for positive po- Figure 3 shows plots of po as a function 
of p for different values of /i and a temperature T — 
120 MeV. As can be seen from this figure, Il(po,p = 0) 
contributes to the increase of the pion mass. Also, for 
large p, the dispersion curves approach the light cone, 
always from within the causal region p\ > p 2 . 

We also notice that the dispersion curves can be 
parametrized by a function of the form 
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(25) 



This is shown in Fig. 4 where the dots represent the com- 
puted values of the dispersion relation and the continuous 
curves the fits. Table 1 shows the values of the parame- 
ters n(T, fj,) and M(T, /x) for the different values of \i and 
the temperature considered. From this table we observe 
that the presence of a finite chemical potential has two 
dramatic effects on the behavior of the pion dispersion 
curve compared to the case where only the temperature 
is considered: first, there is a significant increase in the 
pion mass and second, the index of refraction parameter 
n becomes larger than unity. 

Recall that the magnitude of the pion group velocity 
v g is given as the derivative of the pion energy po with 
respect to the momentum p. From Eq. (Pq) we see that 

dp 
dp 
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therefore, the pion dispersion relation is given as the so- 



where we have defined the (momentum dependent) index 
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fj, (MeV) n M (MeV) 
1 144.411 

120 1.004 150.388 
130 1.006 151.826 



TABLE I: Values of the parameters n(T, and M(T, fi), for 
three values of the chemical potential, fj, — 0, 120, 130 MeV 
and for a temperature T = 120 MeV. 



of refraction n(p) as 

HP) = n , (27) 

and the vacuum pion group velocity w™ c as 




Equation ( f27| ) can be interpreted as a temperature, den- 
sity, and momentum-dependent pion index of refraction. 
Since for finite T and /i, we always have both n > 1 
and po/p™ > 1j f° r au values of p, the index of refrac- 
tion developed by the pion medium at finite density and 
temperature is always larger than unity. 

IV. SUMMARY AND CONCLUSIONS 

In this paper we have considered the effects of a fi- 
nite pion density on the pion dispersion curve at finite 
temperature. The finite density is described in terms of 
a finite pion chemical potential. We have argued that 
such description is valid during the hadronic phase of a 
collision of heavy nuclei at high energies between chemi- 
cal and thermal freeze-out when the strong pion-number 
changing processes have driven the pion number to a 
fixed value. 

In order to consider a general scenario that takes into 
account resummation effects, we have presented an ap- 
proximate solution to the Schwinger-Dyson equation for 
the momentum-dependent pion self-energy, writing this 
as II(P) ~ n +n(P) and considering ft(P) < n , which 
is justified given that in our approximation, II ~ 0(a 2 ), 
whereas the perturbative expansion of IIo starts at 0(pt). 

The pion dispersion relation thus obtained at finite 
density and temperature deviates from the vacuum dis- 
persion relation and can be described in terms of a den- 
sity and temperature dependent mass and an index of 
refraction, also temperature, density as well as momen- 
tum dependent. This index of refraction is larger than 
unity for all values of the momentum. Similar results 
have been obtained in Ref. Q using also a linear sigma 
model without the introduction of a finite chemical po- 
tential but rather as a consequence of the loss of Lorentz 
invariance when a particle travels in a medium at finite 
temperature. We note however that our result is more 



general than that of the former reference where the com- 
putation of the pion dispersion relation was carried from 
the onset in the weak coupling regime, that is to say, 
A 2 C 1, whereas our approach is valid for arbitrary val- 
ues of A 2 @. 

A very interesting consequence of the development of 
such index of refraction happens when considering pions 
that approach the boundary between the strongly inter- 
acting hadronic region of the collision and vacuum. Prop- 
agation through the boundary wall can be described in 
three different regimes, depending on whether the wall's 
width d is much smaller, on the order of, or much larger 
than the pion's wave length A. In the first situation (thin 
wall regime), it is possible to treat the propagation of the 
pions through the boundary in terms of transmission and 
reflection coefficients. In the second case, diffusion effects 
have to be properly taken into account. However, in the 
third case (d <C A), which we choose here as an illus- 
trative scenario, it is possible to consider the motion of 
pions through the wall as that of classical particles climb- 
ing out of a potential well and the interaction of the wall 
on the pions as that of a classical force that causes the 
momentum component normal to the boundary, p t , to 
change, while the momentum component parallel to the 
boundary, pi, is left unchanged 0]. The change is found 
by imposing energy conservation for the waves on both 
sides of he wall 

y/n-i(T,n)tf B + M*(T,vL) = (29) 

which in turn implies 

(p 2 t U = n-\p 2 t ) in - p 2 [l n- 1 } + [M 2 ml] . (30) 

Notice that when n = 1, Eq. ([50|) means that the mo- 
mentum component normal to the boundary in vacuum 
grows when compared to the same component inside the 
pion medium. However, in the general case of finite 
T and /x, < 1 so, from Eq. (|3(]) we see that it 

is possible that for certain values of pi, pf[l — 7i~ 1 ] > 
n _1 (p 2 ) in + [M 2 — m 2 ] and thus that (pt) out be zero or 
purely imaginary, meaning that the pion is reflected back. 

These phenomena might have a word on the behav- 
ior of the transverse pion spectra measured from SIS 
(Heavy Ion Synchrotron) to RHIC energies, but clearly 
a more detailed study of the propagation properties of 
pions through the boundary between the hadronic phase 
and vacuum in relativistic heavy-ion collisions is called 
for. This will be the subject of a future work |20| . 

Acknowledgments 

A. Aranda and A. Ayala wish to thank Universidad de 
Colima for their kind hospitality during the time when 
part of this work was done. Support for this work has 
been received in part by D G APA-UN AM under PAPIIT 
grant number IN108001, CONACyT under grant number 
35792-E and No. 32279-E, the Department of Energy 



7 



under grant DE-FG02-91ER40676 and by Fondo Alvarez- Buylla, Universidad de Colima. 



[1] For a recent review see Quark Matter '01, Proceedings of 

the 15th International Conference on Ultra-Relativistic [13] 
Nucleus-Nucleus Collisions, Nucl. Phys. A698, 1 (2002). 

[2] E. V. Shuryak, Phys. Rev. D 42, 1764 (1990). [14] 

[3] C. Gale and J. Kapusta, Phys. Rev. C 35, 2107 (1987); 
L. C. Liu and W.-H. Ma, ibid 43, R935 (1991). 

[4] J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 [15] 
(1987). 

[5] A. Schenk, Phys. Rev. D 47, 5138 (1993). [16] 
[6] D. Toublan, Phys. Rev. D 56, 5629 (1997). 

[7] M. Loe we and C. Villavi cencio, Thermal pions at finite [17] 

density, [hep-ph/0206294 (2002). 
[8] H. Bebie, P. Gerber, J.L. Goity and H. Leutwyler, Nucl. [18] 

Phys. B 378, 95 (1992). 
[9] P. Braun-Munzinger, J. Stachel, J. Wessels, and N. Xu, 

Phys. Lett. B 344, 43 (1995), ibid 365, 1 (1996). [19] 
[10] CM. Hung and E. Shuryak, Phys. Rev. C 57, 1891 

(1998). [20] 
[11] C. Song and V. Koch, Phys. Rev. C 55, 3026 (1997). 
[12] D. Teaney, Chemical freezeout in heavy ion collisions, 



nucl-th/0204023j (2002). 

L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974); S. 
Weinberg, ibid 9, 3357 (1974). 

A. Ayala, S. Sahu and M. Napsuciale, Phys. Lett. B 479, 
156 (2000); A. Ayala and S. Sahu, Phys. Rev. D 62, 
056007 (2000). 

B. W. Lee, Chiral Dynamics (Gordon and Breach, New 
York, 1972). 

J.I. Kapusta, Finite Temperature Field Theory (Cam- 
bridge University Press, Cambridge, 1989). 
M. Le Bellac, Thermal Field Theory (Cambridge Univer- 
sity Press, Cambridge, 1996). 

H.C.G. Caldas, A.L. Mota and M.C. Nemes, Phys. Rev. 
D 63, 056011 (2001); H.C.G. Caldas, ibid 65, 065005 
(2002). 

R.D. Pisarski and M. Tytgat, Phys. Rev. D 64, R2989 
(1996). 

A. Ayala, work in progress. 



